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Abstract. We study conformal mappings from the unit disk to circular-arc quadrilat- 
erals with four right angles. The problem is reduced to a Sturm-Liouville boundary 
value problem on a real interval, with a nonlinear boundary condition, in which the 
coefficient functions contain the accessory parameters t, A of the mapping problem. 
The parameter A is designed in such a way that for fixed t, it plays the role of an 
eigenvalue of the Sturm-Liouville problem. Further, for each t a particular solution 
(an elliptic integral) is known a priori, as well as its corresponding spectral param- 
eter A. This leads to insight into the dependence of the image quadrilateral on the 
parameters, and permits application of a recently developed spectral parameter power 
series (SPPS) method for numerical solution. Rate of convergence, accuracy, and com- 
putational complexity are presented for the resulting numerical procedure, which in 
simplicity and efficiency compares favorably with previously known methods for this 
type of problem. 
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Introduction 

A symmetric circular quadrilateral (s.c.q.) is a Jordan curve P in C formed of 
four circular arcs (or straight segments) with all four internal angles equal to 
7r/2. Will assume that the vertices of P are situated in positions of the form 
±A, ±A, where Re A > 0, Im A > 0. 

Let D be a plane domain containing the origin and bounded by an s.c.q. The 
set of conformal types of such domains (considering the vertices as distinguished 
points) forms a two-dimensional space in a natural way. We consider conformal 
mappings /:D — > D from the unit disk D to D. Since / extends continuously 
to the boundary P, there is a unique value t e [0, 27r] such that /(e**) = A; We 
will generally assume that /(O) = and /'(O) > 0, which implies < t < n/2. 

The more general problem of mapping the disk to circular-arc polygon do- 
mains is treated in [^, [DTl Chapter 4], [Hil Chapter 16], ^ Section 17.6], 



1 Partially supported by CONACyT grant 50424 

2 Partially supported by CONACyT grant 80503 



1 



|Ho| . |Ne[ Chapter 5]. In particular, it is well known that the Schwarzian deriva- 
tive of a conformal mapping of D onto a circular-arc polygon is a rational func- 
tion of degree two. In this article we will develop further the work of P. Brown 
on the accessory-parameter problem for s.c.q.s. We follow much of the notation 
and copy several equations from [B] where it is verified that due to the symme- 
tries of s.c.q.s, the Schwarzian derivative Sf of / is of the specific form ([1]) below. 
This function is determined by two real parameters s and the relationship ([2]) 
must hold. As a partial converse, it is well known that due to the intimate re- 
lationship between the Schwarzian derivative and curvature, if the Schwarzian 
derivative of a holomorphic function / defined in D is of the form (P) and if © 
holds, then / is a local homeomorphism onto a (not necessarily schlicht) domain 
bounded by a (not necessarily simple) right circular-arc polygon. 

A basic question is the following. Given P — dD (for example by specifying 
A and also the radius or the midpoint of one of the edges of P), to find the 
parameters t, s of the Schwarzian derivative of /. In [B] Brown looked first at a 
simpler question, fixing t, normalizing /'(O) = and calculating the remaining 
parameter s corresponding to a geometric characteristic of P, such as the radius 
of one of its edges. This reduces the problem to one real dimension. As with 
most methods which have been developed for conformal mapping of circular-arc 
polygons, in [B] the Schwarzian differential equation is solved numerically for 
trial values of s, the corresponding geometric characteristics of P = /(OD) are 
calculated, and the process is repeated until a sufficiently close s is located. 
Here we apply a technique developed in |KPj for dealing with Sturm-Liouville 
problems, the spectral parameter power series (SPPS) method, which permits 
a direct calculation in the sense that once certain auxiliary parameters are cal- 
culated for given t, one may evaluate the solution corresponding to any desired 
s without resorting to further integrations. 

In a second article |B]] the full two-parameter problem is addressed. The 
problem is formulated with the normalization /(I) = 1, and the data is given in 
terms of the radius 1/ki of the right edge of P and the midpoint P2 of the upper 
edge. Brown's solution involves a table of previously calculated values of these 
parameters in terms of {t, s). To apply it one looks for values reasonably close 
to the desired geometric parameters in this table, and then applies an iterative 
process to approximate the sought-after {t, s) to the desired accuracy. In this 
paper we apply our solution of the one-parameter problem, which is quite rapid, 
to this two-parameter problem in an iterative way. In part due to properties 
of an equivalent parameter A which we use in place of s, our method does not 
require consultation of a table of prior values. 

In the next two sections we set up a Sturm-Liouville boundary value problem 
whose solution relates the accessory parameter A to the curvature k of the right 
edge of P. It is seen that A is a spectral parameter in a boundary value problem. 
In Section[3]we describe the so-called canonical mapping f^o of D to a rectangle, 
and identify its parameter Aoo and the corresponding eigcnfunction ?/oo of the 
Sturm-Liouville problem, which are used in the application of the SPPS method 
which is summarized in Section |4] and then applied in Section [5] to represent k 
as a power series in A for fixed t. An algorithm for the one-parameter problem 
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K X for fixed t is presented in Section [SI and for the two-parameter problem 
(ki,P2) in Section [71 



1 Sturm-Liouville equation 

We begin by setting up the classical second-order differential equation which 
governs the conformal mapping to s.c.q.s. The Schwarzian derivative of a holo- 
morphic /, namely 5/ — {f" /f')' — (1/2) (/"//')^, is again holomorphic when 
/' does not vanish. For the conformal mapping of D to an s.c.q., formula [BJ 
(3)], which we will refer to as (B3), says 

e / ^ 3 / 1 1 1 1 



c c c c 



z + e 



(1) 



in which the parameter c G C is subject to one restraint as follows. Write 
c = pe''. Then (B8) 

_3 

^ ~ 8 cos(s -f- 1) ^^-^ 

where further we may assume (B9) Tr/2 — t < s < 3tt/2 — t. 

We now introduce a new parameter equivalent to s. By the quantity 
p > satisfies in fact p > 3/8. Since s = arccos(— 3/(8/9)) — t where tt/2 < 
arccos(— 3/(8/?)) < 37r/2, the parameter c in (|T]) can be expressed as follows: 

c = pe^'' = p ( cos((arccos — ) — t) + i sin((arccos — ) — t) 
V 8p 8p 



= — {cost — i sin t) + e^yp^--(3JW)^ {sint + i cost) 
8 
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--e 



-it I \-„-it 



Aze"** (3) 



where we write 
with 



= e^p^ - (3/8)2 

-1, s < n -t, 

0, s = IT — t, 

1, TT — t < S. 

At the particular value s = tt — t we have p = oo, and also p — > cx) as s tends to 
either extreme value tt/2 ~ t or 37r/2 — Further, it is easily seen that 

3 

A=-tan(s + t). (4) 
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Everything we will do depends on the fact, manifested in ([3]), that the pa- 
rameter c is a linear polynomial in A. The same then holds for Sf{z). Indeed, 
abbreviating a = e**, we find that 

o / N / 1 1 1 



--( — 

4 V z2 - 



1 



+2Xi 



1 1 



where 



2 V(22-a2)2 (22_^2)2^ 
^l(^) = ~M .2 ^2 ~ ^2 ^ ^2 ) (5) 



From now on we will regard Sf as parametrized by t, A instead of t, s. The 
functions ipQ^ ipi depend also on t, so when we require more clarity we will write 

Sf{z) — Rt^\{z) where 

i?t,A = 2Vo-2A0i. (6) 



It is a textbook fact |Hil INe| that all holomorphic functions having Schwarzian 
derivative equal to Rt,\ are quotients 2/2/2/1 where j/i, 2/2 are linearly independent 
solutions of the ordinary differential equation 2y"+Rt.\ i/ = in ©, i.e., solutions 
of 

y" + -ipoy = A?/'i2/. (7) 

It is our purpose to take advantage of the particularly nice form of equation ([7| 
and its relation to spectral theory. 



2 Boundary Value Problem 

In [B] Brown describes in detail the possible relations among the edges of s.c.q.s. 
In particular, for a generic s.c.q., the two circles containing one pair of opposite 
edges will intersect, while the two circles containing the remaining two edges 
are disjoint. There are two degenerate cases: when one pair of opposite edges 
lie on the same circle, or when all four edges are straight segments (i.e., when 
P is rectangular). For the generic case, following [B we will take the right and 
left edges to lie in disjoint circles of common radius |r|, where r is given (B14) 

by 



/'(I) + /"(!)■ 



(8) 
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We will normalize / by 



/(0) = 0, /'(0)-l, /"(0)-0 



(9) 



recalling that the Schwarzian equation has three degrees of freedom. (This 
means that the image may need to be rescaled, i.e., P — constant • /(9D)). 

Let y be a solution of ([7]) in D normalized by y{Q) — 1, y'(0) ~ 0. Then the 
normalized mapping / is equal to the indefinite integral 



/(^) = 



yicy 



■dC. 



(10) 



Indeed, ((TO)) yields f — y and /" — —2y ^y', so / satisfies the required 
normalization (l9]). Now from (El) we have 



= r 



( /'(I) + ,/"(!) ) - /'(If = r ( yil)-' - 2y(l)- V(l) ) " y(l)- 



so the curvature of the right edge of the image is given by 

« = -=2/(l)'-2y(l)y'(l). 



(11) 



Consequently the boundary value problem to solve is for y{z) with real z e [0, 1], 

'' y" + V'oy = AV'iy, 

2/(0) = 1, 
2/'(0) = 0, 
L y{lf -2y{l)y\l) = K. 



(12) 



The data for this problem are t and k, and one seeks A for which there is 
a solution to the system (fT^ . where we recall that the dependence on t is 
through Tpo a-iid ipi . Among such A one of course wants the value for which the 
corresponding mapping / is injective. In ^ there are many details regarding 
the values of t, s for which this holds, and one can apply these facts to A. 

The nonlinear nature of the boundary condition at z = 1 in (jl2p will be 
dealt with in Section O 



3 The Degenerate S.C.Q.s 

Again following [B], let 0{s) S M denote the center of the circle containing the 
rightmost edge of the s.c.q. /(D). In [B^ it is shown that 0{s) is a monotone 
function of s for each fixed t. Two special values sq, Sqo are singled out, with 
O(so) — (called the root mapping) and O(soo) — oo (the canonical mapping), 
for which /(D) has the degenerations mentioned at the begining of Section [51 
The root mapping and its parameter sq are worked out in detail. We will use 
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Soo here because it is simpler; in fact, the canonical mapping fao is an elliptic 
integral of the first kind: 

foo{z)= r(l-2(cos20C' + C')-'/'dC (13) 
Jo 

Here /oo(B') is a rectangle, /oo([— 1,1]) ^ K, and the values ki — /oo(l) and 
^2 — —ifoo{i) are known clasically as complete elliptic integrals. 

Note that f'^ never vanishes, and /^(O) = 1. By direct calculation one finds 
that 



Sf 



f'oo) 2 

_ 2(cos 2t + (cos^ 2t - ?,)z^ + (cos 21)7^) 
(1 - 2(cos2i)^2 + ^4)2 • 

However, by (B23), 

cos(s — i) 

_ 2(cos2i + (cos2 2t- 3)z2 + (cos2t)z'') cos(s + t)^ 

~ (l-2(cos2i)z2 + z4)2 + 2(l-2(cos2t)z2 + z4)- ^ ^ 

By comparison, the second term on the right side of (|14p must vanish, 

cos{s-t) 
3 — cos2t. 

cos(s + t) 

This can be written many ways, for example 

,cos2t-3 , ^-rir+sint 

tan s = (cot t) = - (cot t) (15) 

^ Vos2t + 3 ^ ^ 1 +cosi ^ ^ 

cos t 

etc. Given t, (|15p provides s = Soo for the canonical mapping. With this in 
hand, we obtain the parameters 

-3 



8cos(soo + 1) ' 



-oo\ pIo 



3^2 



-3 

— tan(Soo + tj- 



With the aid of (|T5)) one can simplify this last equation to 



Aoo = -cot2i (16) 
4 ^ ' 
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The Sturm-Liouville system (jl2p corresponding to the canonical mapping 
with Schwarzian derivative Rt,x^ has the particular solution j/oo = (/oo)~^^^) 
specifically 

yoo(z) = (l-2(cOS2t)z2 + z4)i 

Clearly this satisfies the normalizations yoo(O) = 1 and y^(0) = 0, and for 
reference we note 

y„,(l) = V2^t, y'^il) ^ (18) 

Since /oo(D) is a rectangle, the curvature of its right edge is — 0. Further, 
since f'^ never vanishes, it follows that the function j/oo also never vanishes in 
D; this can also be seen from (|T7| . 



4 Solution by Iterated Integrals 

We now apply the spectral parameter power series (SPPS) method for solution 
of Sturm-Liouville problems, developed in [KPj . This will provide power series 
in A which represent the geometric parameters of the s.c.q. mapping problem. 

Let two functions qo, qi be given on [0, 1]. (All of the following is valid for 
qn{z) defined in, say, \z\ < 1, but the numerical method presented later will not 
involve complex values of z.) The sequence /„ of iterated integrals determined 
by the generating pair ((?0i9i) is defined as follows. Let Ii = 1 identically on 
[0, 1], and then recursively for n >2, 

In{z)^ I In-liOqn-liOdC (19) 

Jo 

where the indices in the generating pair Rre understood mod 2; i.e., Qn-\-2j ^ Qn 
for j = 1,2,... 

The relationship between the iterated integrals and the basic equation ([7]) is 
given by the following result. 

Proposition 4.1 |KP| Let ipo and ipi be given, and suppose that yaa is a func- 
tion which does not vanish and which satisfies the ordinary differential equation 

on [0,1]. Choose qo = 1/j/^, qi — ipivlo '^'^'^ define X^"), X^") to be the 
iterated integrals determined by {qo,qi) and by {qi,qo) respectively. Then for 
every A € C the functions 

oo 

yi^yooY.i\-Xoofx(''\ 

oo 
k=0 
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are linearly independent solutions of the equation 



on the interval [0, 1]. Further, the series for yi and 7/2 converge uniformly on 
[0, 1] for every X. 

5 Variation of Curvature 

To apply Proposition 14.11 to the boundary value problem we take ipo,i^i 
as in ([5]) and yoo as in (fT7|) . Note that 



with a similar formula for y'2 , so the initial values of the particular solutions yi , 
2/2 are given by 



This says that yi as given by Proposition 14.11 is precisely the solution of the 
second order linear differential equation which satisfies the two initial conditions 
of at the point z — 0. Consequently problem reduces to satisfying the 
final boundary condition, 




J/i(0) = yoo(0) = l, y[{0) 
2/2(0) -0, 2/^(0) 



yoo(o) 



2/i(l)(2/i(l)- 22/1(1)) 



(20) 
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-25 



Figure 1: Graph of function of A for f = O.Itt, 0.2tt, O.Stt, OAtt. 



8 



The numerical evaluation of the left side of ([^0]) is readily accessible because it 
is a power series an{X — Aoo)" in A whose coefficients a„ are represented in 
terms of those of the series yi(l) = J2 ^n(A— Aoo)" provided by Proposition l4.1l 

In Figure [T] we show a graph of the power series (|^D|) for several values 
of t. Numerical details will be given later; here we wish to stress that once 
having calculated the coefficients a„ for a given t, one needs no longer to solve 
a (Schwarzian or Sturm-Liouville) differential equation in order to calculate k 
as a function of A for that value of t. 

Perhaps surprisingly, the calculation of (|20p can be simplified much further. 
Note that the k-th Taylor coefficient of yi{l) — 2y'i{l) is 

bk~2b', = 2;^(l)x(2'=)(l)-2fy:^(l)X(2'=)(l) + ^X(2'--i)(l)) 

V 2/00(1) J 

2/00(1) 

It follows from ^T7\ that ?/oo(l) ~ 2y^(l) — 0, which leaves 
bk-2b', = L^X(2'^-1)(1). 

yoo(i) 

Thus 

n 

a„ = -2 ^ X^^/^) (1) (1). (21) 

From ((5|l. (fT7ll it is also clear that yoo{z) > and ipi{z) > for z G [0, 1], and 
consequently X^'^^z) > 0. Thus we see that in fact all the Taylor coefficients 
a„ of K as a function of A are negative, with the exception of aq = 0. 

Note that k is never zero for A > Aoo! m fact, dn/dX is strictly negative, 
and ooasA^oo. Thus the radius of the rightmost edge of /(D) tends 

to zero, as illustrated in Figure [H As A — > 00, the arc e^^ : —t < 9 < t oi 




Figure 2: Images /(D) for t = 7r/4 and A — 1.3, 1.4, 1.5. The images are depicted 
near the point w = 0.5, where it is perceived that the rightmost figure is a non-schficht 
region. 
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-1.5 t 



Figure 3: Images /(D) for t = 7r/4 and (a) A = = A^o = 0, (b) A = -0.32219, 
A — —0.91570 and (d) A = —1.43554 (schematic drawing of boundary). The value in 
(d) corresponds to second largest zero of fi:(A), after Aoo. The values for (b), (c) were 
chosen to give equal values k = 0.8. 

f{dD) is mapped to a curve tracing ever-greater numbers of full turns around 
ever-smaller circles. 

In Figure [3] we bend the edge in the opposite direction, taking A < Aoo- The 
largest zero of k(A) is Aoo corresponding to the canonical mapping. The next 
zero gives a domain, exterior to the non-simple curve depicted in Figure [3fd), 
covering the exterior in the Riemann sphere of two circles joined by a straight 
segment. Further zeros correspond to non-schlicht domains. 

6 Algorithm for the One-Parameter Problem 

Here we solve the nonlinear problem mentioned at the end of Section [2l 
Fix < t < tt/2. Given k, we seek a parameter A so that the mapping / with 
Schwarzian derivative Sf = Rt.x, normalized by /(O) = 0, /'(O) = 1 will send 
dD to an s.c.q. whose right edge has curvature k. (If k < 0, we want the center 
of the edge arc to lie to the left, i.e., in the ray (— oo, 1).) We stress again that 
for fixed once certain parameters have been determined, the relation k i— > A 
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is of rapid calculation. 

The algorithm given below summarizes results of the previous sections, in 
particular formulas (O, (fT7|) . and (I^T]) . In the following, "evaluate" means to 
calculate at M + 1 evenly spaced values of z from to 1 . We will use iV + 1 
terms in the power series in A. Observe that V'o a-nd X*^"^ are not needed. 

One-parameter algorithm. 

1. Evaluate yoo{z) = (1 - 2(cos2t)z2 + ^4)1/4^ 

2. Evaluate 

2sin2t 

"^'^^^ " z4-2(cos2t)z2 + l- 

3. Evaluate the iterated integrals X^"^ determined by (V'ly^, for 
n = 0,l,...,2N. 

4. Calculate the Cauchy product terms 

n 

a„ = -2^1(2fc)(^)^(2(„-fc)-i)(^^_ 

for n = 0, 1, . . . , TV. Calculate also Aoo = (1/4) cot2t. 

5. For any desired k, solve the polynomial equation 

N 

^ a„(A - Aoo)" = K 

n=0 

for A — Aoo, and then add Aoo to the result obtain A. 

Then the mapping ft^x will produce curvature approximately k on the right and 
left edges. 

7 Two-Parameter Problem 

We now turn to the two-parameter mapping problem addressed in [B] . Given an 
s.c.q. P proposed to be the image of D under w — f{z), one does not know the 
value of /'(O) directly from the geometry of P, so instead one assumes that the 
right edge of P passes through w = 1. Consider the mapping g{z) = /(z)//(l) 
which satisfies 

5(0) = 0, 5(1) = 1- 

We let Ki and K2 denote the curvatures of the right and upper edges of P = 
g{dlD)), with midpoints pi = 1 and p2 G iM^. The basic mapping problem will 
be to find (i. A) given (ki,P2)- 
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Geometric parameters of P. To begin with, we take into account the values 
which we could readily calculate if we were already given {t, A). Write wi — /(I), 
W2 = /(«)• We have already seen how to calculate wi = 2/2(l)/2/i(l)- For W2, 
observe that by ([51), ([11), 

K/2-t,^x{±^z) = -RtA^)- (22) 

Further, u'2/j — where f*{z) — —ifliz), and Chain Rule for Schwarzian 

derivatives gives us 

Sf{z) = -Sf{iz) = Rv,\'{z). 

Thus W2 can be obtained simply by using (7r/2 — i, —A) in place of (i. A) in the 
calculation of wi . This provides us in turn with the desired values 

W2 

P2 = 

Wi 

and 

Kl = KWl 

where k denotes as previously the curvature of the right side of {l/wi)P = /(D)- 
As we have seen, for fixed t these quantities are power series in A. 

Calculation based on t and ki . Now we consider how to calculate A correspond- 
ing to given values of t and ni . We know that 

Ki = ^/(1) = ^^, 

yi(i) 

A graph of this function of A is illustrated in Figure ID Since the right edge 
of /(HJ) passes through 1, the radius of this edge can never be less than 1/2. 




-140 -120 -100 -80 -60 -40 -20 



Figure 4: Normalized curvature ki as a function of A for t — O.Stt 
through the local maxima at height ki = 2, the boundary image 
behavior shown in Figure [3jd). 



As A passes 
exhibits the 
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Thus we need only look for values of A greater than the largest value for which 
Kl = 2. 

An obvious method of calculation is to try different values of A, using (i, A) 
to calculate ki , and by bisection or another similar method to close in on the de- 
sired value of Ki. However, one can take advantage of the equivalent formulation 
a(A) = where 

a{\)^nv2{l)~Kiv^{l). (23) 

The coefRcients of the power series a(A) are found in terms of known series k, 
2/1, and t/2 in A (note that the first term on the right side of (^5)1 is a Cauchy 
product). Of course, the calculation of 2/2(1) involves the iterated integrals 
which were not used in the 1-parameter algorithm. Thus the numerical 
calculation of A reduces to finding the zero of a polynomial with least absolute 
value. 

Finally we consider the upper midpoint p2- 

Lemma 7.1 For any fixed value of ki G M, the geometric parameter p2 is an 
increasing function of t ^ (0, 7r/2). 

Proof. Fix K2, and let t < t' . The topological quadrilateral bounded by 
the unit circle dD with vertices ±e^'* has conformal module smaller than the 
one with vertices ±e^** . Likewise, the s.c.q. P whose upper edge meets the 
imaginary axis in iy has conformal module smaller than the one meeting in iy' , 
when y < y' (see Figure [5]). The statement follows. □ 

Two-parameter algorithm. We calculate (i. A) from (ki,P2)- Choose t~ and t~^ 
near to and 7r/2, respectively, so that the desired t can be sought for in the 
interval < t < i+. Begin with the midpoint t — {t^ +i+)/2, and for the three 
values t~ , t, t^ each paired with ki, find the three corresponding values of p2 
as described above, say p, p^ . By Lemma [SJ we may repeat the process in 
whichever of the intervals (p^,p) or (p,p^) the desired p2 may be located. 



iy 





n/ 








/fAe"') 
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Figure 5: Increasing conformal module of P corresponding to increasing p2. 
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Figure 6: Images /(D) exhibiting ki = —1 and p2 — 1.0, 2.0, and 3.0 respectively. 

In Figure [6] some examples are given of s.c.q.s with prescribed values of ki 
and p2, calculated via this two-parameter algorithm. 

Variations on the above considerations can easily devised for finding (t, A) 
in terms of (ki,K2), etc. 

8 Domain of univalence. 

The above considerations permit us to calculate easily the complete set of (t, A) 
for which the solution to Sf = Rt^\ is univalent on ©. From the discussion 
of Figure [21 as A decreases starting from Aoo, the uppermost vertices /(e**), 
/(e'f'^/^"*)) meet on the positive imaginary axis at a certain critical value 
A = Ajiiin. At this value the right and left image edges, being orthogonal to 
the upper edge, must lie within the extended imaginary axis. From this it fol- 
lows that /(I) = oo, /(—I) — — oo, and consequently yi{l) ~ 0. (Note that 
2/2(1) by hnear independence.) Of course by (^0)1 this implies k = 0, as 
was previously discussed. Similarly, as A increases from Aoo the critical value 
Amax approximated in Figure [2] is attained when /(e**), /(e"**) meet on the 



A 




4 2 



Figure 7: (left) Domain of values of {t, A) for which the conformal mapping to an 
s.c.q. is univalent, (right) Domain rescaled by application of arccot 4A to the vertical 
axis. 
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positive real axis, and f{i) = ioo, f{—i) = —ioo, so instead we solve yi{i) — 0. 
In view of ((22)) this can be solved by using 7r/2 — t in place of t. Thus for each 
t e (0,7r/2) we have values 

Amin 

(t) < Aoo(i) < A,„ax(i) 

and in fact 

Amin(7''/2 — t) = — Ainax(i)- 

The domain of univalence {{t, A): Xmin{t) < i < Aniax(0} shown in Figure [7l 

9 Numerical Results 

Calculations were performed in Mathematica 6.0 (Wolfram). 

The numerical integrations were carried out by subdividing [0, 1] into groups 
of 5 consecutive intervals of length 1/M, and multiplying each group by a ma- 
trix equivalent to integrating the degree-5 polynomial passing through the cor- 
responding values of the integrand. 

In the calculation of the iterated integrals /„ determined by functions 50, 
qi satisfying a bound |(7„(z)| < from the definition we have |/„(z)| < 
i^(sup |/„_i|)2; for z e [0,1], and thus \In{z)\ < K^z^/n\. From this it is 
seen that the Taylor coefficients of the power series we are considering tend to 
zero quite rapidly. Instead of working out here the simple a priori error bounds 
which result from this, we follow the time-honored procedure of observing how 
the numerical method actually works. In terms of the two main parameters, 
M (the number of subdivisions of [0, 1] and N (the degree of the polynomial 
approximation to the power series), the following tables show the precision ob- 
tained for different values of {t, A). Naturally, as | A — Aoo| increases, so does the 
numerical error (recall that each Aqo tabulated depends on the corresponding 
value of t). 

Required values of M, N for calculation of k to 5 significant figures 



A== 


Aoo " 


- 1 


A = 


Aoo - 


- 2 


A = 


Aoo ~ 


- 5 


t/TT 


M 


TV 


t/lT 


M 


iV 


t/lT 


M 


N 


0.05 


40 


15 


0.05 


40 


15 


0.05 


60 


15 


0.08 


25 


15 


0.08 


25 


15 


0.08 


25 


15 


0.1 


15 


15 


0.1 


20 


15 


0.1 


20 


15 


0.2 


10 


15 


0.2 


15 


15 


0.2 


20 


15 


0.3 


10 


15 


0.3 


15 


10 


0.3 


20 


15 


0.4 


10 


11 


0.4 


10 


11 


0.4 


15 


15 


0.45 


10 


11 


0.45 


10 


11 


0.45 


15 


15 


0.48 


10 


9 


0.48 


10 


9 


0.48 


10 


12 



15 



Required values of M, N for calculation of k to 8 significant figures 



A = 




- 1 


A = 




- 2 


A = 




5 


t/n 


M 


N 


t/n 


M 


N 


t/n 


M 


N 


0.05 


85 


25 


0.05 


90 


25 


0.05 


110 


35 


0.08 


75 


25 


0.08 


65 


25 


0.08 


110 


35 


0.1 


65 


25 


0.1 


50 


25 


0.1 


110 


35 


0.2 


40 


25 


0.2 


30 


25 


0.2 


85 


35 


0.3 


30 


20 


0.3 


25 


35 


0.3 


40 


30 


0.4 


25 


15 


0.4 


25 


20 


0.4 


30 


25 


0.45 


25 


15 


0.45 


25 


15 


0.45 


25 


20 


0.48 


20 


15 


0.48 


25 


5 


0.48 


20 


15 



One surprising fact is that the algorithm has no difficulty approaching t = 
7r/2, even though the conformal module of the s.c.q. tends to infinity and the 
well-known crowding phenomenon makes graphing difficult. 

In the 1-parameter algorithm, steps 1 and 2 have computational complexity 
of order 0(M); step 3 is 0{MN) and step 4 is 0{N'^). The solution of roots of 
polynomials with real coefficients is now a well-refined science which we will not 
expound upon here; one can use with goods results the internal function NRoots 
of Mathematica, which we believe has an operation count of the order 0{N'^). 
However, it is even simpler to use the standard routine FindRoot for finding 
a single zero of an arbitrary function, giving Aoo as the starting point of the 
search. Each evaluation of the function requires 0{N) arithmetic operations, 
and the number of evaluations (by a modified form of Newton's method) is less 
than proportional to the logarithm of the accuracy sought. 

While the parameters t, A have been obtained by the SPPS method, the cal- 
culation of the image domains in Figures [2l [3l and [6] are based on the function 
NDSolve, which was used to integrate the second-order linear differential equa- 
tion along radii of ©. This is solely for illustration of the image domains, not for 
application of the algorithms presented here. In these figures the vertices, being 
singular images, were determined in terms of neighboring boundary points. 

The two-parameter algorithm gives no essentially new numerical considera- 
tions. The number of iterations of the bisection process corresponds exactly to 
the desired accuracy, which also determines the choice of M, iV as reflected in 
the tables above for the calculation of k. 

10 Conclusions 

p. Brown (BJ [B] studied the accessory parameters which govern conformal map- 
ping from the unit disk to a symmetric circular quadrilateral, and gave algo- 
rithms for finding these parameters in terms of geometric characteristics of the 
quadrilateral. 

We have introduced an alternative parameter A in the expression for the 
Schwarzian derivative of the conformal mapping, and shown that it is identified 
naturally with a spectral value of a Sturm-Liouville problem with a nonlinear 
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boundary condition. The SPPS method in |KP| for solving Sturm-Liouville 
differential equations provides simpler algorithms, especially when the solution 
for a particular value of the spectral parameter is known a priori, as is the case 
for the problem considered here. 

One advantage of the SPPS approach as compared to other methods is that 
once the simple procedure of calculating the indefinite integrals is carried out, 
the problem reduces to finding roots of a polynomial. To find the accessory 
parameters it is not necessary to evaluate the conformal mapping explicitly at 
interior points of the domain. Further, the method presented here does not 
require locating the desired geometric characteristics approximately in a table 
of previously calculated values as was done in [B] . 

Because of the natural relation of the second-order ordinary differential equa- 
tion to the Schwarzian derivative, we believe that the SPPS approach given here 
may be applied in a similar way to a wide variety of conformal mapping prob- 
lems. 
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